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Abstract 

This paper presents a sampled-data approach for the robust control of a single qubit (quantum 
bit). The required robustness is defined using a sliding mode domain and the control law is designed 
offline and then utilized online with a single qubit having bounded uncertainties. Two classes of 
uncertainties are considered involving the system Hamiltonian and the coupling strength of the system- 
environment interaction. Four cases are analyzed in detail including without decoherence, with amplitude 
damping decoherence, phase damping decoherence and depolarizing decoherence. Sampling periods are 
specifically designed for these cases to guarantee the required robustness. Two sufficient conditions are 
presented for guiding the design of unitary control for the cases without decoherence and with amplitude 
damping decoherence. The proposed approach has potential applications in quantum error-correction and 
in constructing robust quantum gates. 
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Nomenclature 

Throughout this paper, we use the following notation: 

| \\f) state vector (quantum pure state) 

a* complex conjugate of a 

L T transpose of L 

l) adjoint of L 

tr(A) trace of A 

{\ff\ adjoint of |y/) 

(0| i//} inner product of \<j>) and \\]f) 

p density operator 

<J x , y ,z Pauli matrices 

CO(t) uncertainty amplitude in o z 

£ x (t) uncertainty amplitude in o x 

£y(t) uncertainty amplitude in a y 

R set of real numbers 

7 coupling strength 

8y t uncertainty in coupling strength 

& c sliding mode domain of closed systems 

3> a sliding mode domain of quantum systems with amplitude damping decoherence 

3> p sliding mode domain of quantum systems with phase damping decoherence 

^ sliding mode domain of quantum systems with depolarizing decoherence 

C coherence 

P purity 

po probability of failure 



I. INTRODUCTION 

Controlling quantum phenomena is becoming an important task in different research areas 
such as quantum optics, physical chemistry and quantum information [Q]|-[|4). The development 
of quantum control theory can provide systematic methods and a theoretical framework for 
analyzing and synthesizing quantum control problems. Several theoretical tools and design 



methods in classical control have been applied to the quantum domain. For example, Lie groups 
and Lie algebras have been used to establish controllability conditions for closed quantum 
systems Q. Optimal control theory has been applied to control analysis and the design of several 
quantum control tasks such as population transfer with minimum energy or in the shortest time 
10- [HI. Learning control has become a powerful tool for the direct laboratory discovery of 
laser pulses controlling a variety of atomic and molecular phenomena 0. Feedback control has 
been utilized for the control of quantum entanglement, quantum error-correction and quantum 
state preparation ffT0l - fl20l . The development of quantum control theory needs to consider the 
special characteristics of quantum systems (e.g., measurement collapse and non-commutative 
relationships) and the unique objectives of quantum control (e.g., entanglement generation and 
decoherence control) (For more discussion, see, e.g., [JT]). 

Robust control is one of the most important research areas in classical control theory. Attaining 
robust control for quantum systems has been recognized as a key issue in the development 
of practical quantum technology [1211 - 11251 . since many types of uncertainties are unavoidable 
(including control noise, environmental disturbances, etc.) for most practical quantum systems. 
Several methods have been proposed for the robust control of quantum systems. For example, 
James et al. If26l formulated and solved a quantum robust control problem using the H°° method 
for linear quantum stochastic systems. A risk-sensitive control problem has been solved for a 
sampled-data feedback model of quantum systems ll27l . Quantum robust control is still in its 
infancy, and it is necessary to develop new tools to deal with different types of uncertainties. 

Dong and Petersen Il28~l - ll30l developed sliding mode control to enhance the robustness of 
quantum systems. In particular, two approaches based on sliding mode design OTI have been 
proposed for the control of quantum systems, and potential applications of sliding mode control 
to quantum information processing have been presented [28] . Sliding mode control for two-level 
quantum systems was presented to deal with bounded uncertainties in the system Hamiltonian 
|]29l . This paper will employ the concept of a sliding mode domain to define the required 
robustness and develop a new sampled-data design approach [|32l . Il33l to enhance the perfor- 
mance of a controlled quantum system with uncertainties in the Hamiltonian as well as in the 
system-environment interaction. 

Sampled-data control has been widely applied in industrial electronics, process control and 
signal processing [|33l . The sampled data are used to design controllers while the sampling 



(measurement) process is usually assumed not to affect the system's state. However, in quan- 
tum control, the sampling process unavoidably destroys the system's state according to the 
measurement collapse postulate (see, e.g., 01). Hence, measurement can be used as the means 
for information acquisition as well as a control tool. For example, several incoherent control 
schemes have been presented where measurements are used as a control tool to affect the system 
dynamics ll34ll - 11361 . A framework of quantum operations including unitary control and projective 
measurements has been developed to investigate feedback control of quantum systems [1371 . Il38l . 
One well known example where measurement modifies the system dynamics is the quantum Zeno 
effect, which is the inhibition of transitions between quantum states by frequent measurement 
of the state (see, e.g., [|39ll and Il40l0 . However, it is usually a difficult task to make frequent 
measurements with practical quantum systems. We may assume that the smaller the measurement 
period is, then the bigger the cost of accomplishing the periodic measurements becomes. Hence, 
in contrast to the quantum Zeno effect, in this paper we will use the sampling (projective 
measurement) process as a control tool and design sampling periods as large as possible to 
guarantee the required robustness for several classes of quantum control tasks including control 
design for quantum systems with uncertainties in the system Hamiltonian and robust decoherence 
control of Markovian open quantum systems. 

Decoherence occurs when a quantum system interacts with an uncontrollable environment BP . 
Decoherence has been recognized as a bottleneck for the development of practical quantum 
information technology [|42l . Various methods have been proposed for decoherence control 
including quantum error-avoiding codes B3l - [|45l . quantum error-correction codes [|46l. dy- 
namical decoupling HTll . Il48l and quantum feedback control ll49ll . In quantum error- avoiding 
codes, quantum information is encoded in a decoherence free subspace which is inherently 
immune to decoherence due to specific symmetries in the system-environment interaction [|50l . 
Quantum error-correction codes are active methods to detect and counteract the effects of errors 
during quantum information processing via encoding redundant qubits. Dynamical decoupling 
of decoherence control is an open-loop control approach which often employs bang-bang control 
pulses to dynamically cancel the effect of decoherence. Quantum feedback and optimal control 
theory also provide powerful tools for the analysis and design of decoherence control [|5D . [|52l . 
However, there are few results which consider robustness when uncertainties or inaccurate param- 
eters exist in the system Hamiltonian or the system-environment interaction. Here we consider a 



robust decoherence control scheme for quantum systems subject to Markovian decoherence [HJ. 
In particular, we will focus on a single qubit subject to amplitude damping decoherence, phase 
damping decoherence and depolarizing decoherence [[4). We propose a sampling -based design 
approach to guarantee the robustness of a single qubit system with uncertainties in the system 
Hamiltonian and the coupling strength of the system-environment interaction. 

The paper is organized as follows. Section [I] presents the control problem formulation and 
defines the required robustness. In Section |nil we present the main methods and results for 
robust control design. Section [IV] gives the proofs of the main results. Concluding remarks are 
given in Section [V] 

II. Control problem formulation 

For an open quantum system, its state is described by the positive Hermitian density matrix 
(or density operator) p satisfying trp = 1, and the evolution of p cannot generally be described 
in terms of a unitary transformation. In many situations, a quantum master equation for p(t) 
(or p t ) is a suitable way to describe the dynamics of an open quantum system. One of the 
simplest cases is when a Markovian approximation can be applied under the assumption of a 
short environmental correlation time permitting the neglect of memory effects BT1 . For an N- 
dimensional open quantum system with Markovian dynamics, its state p(t) can be described by 
the following Markovian master equation (for details, see, e.g., iPffl . [[531 , 11541 ): 

p(o = -/[//(o,p(0]+^E 1 a^[^pw,4]+[ L i 5 p(o4]}- (d 

j,k=0 

Here for an arbitrary operator X, [X,p] =Xp — pX is the commutation operator, {L/}y = Q is a 
basis for the space of linear bounded operators on the Hilbert space with Lq = I, the coefficient 
matrix A = (ot,^) is positive semidefinite and physically specifies the relevant relaxation rates and 
we have set % = 1 in this paper. Markovian master equations have been widely used to model 
controlled quantum systems in quantum control [I55l - [l57l . especially for Markovian quantum 
feedback 0. 

In this paper, we will focus on a two-level quantum system (a single qubit) with Markovian 
dynamics whose evolution can be described by the following Lindblad equation: 

K 

p(t) = -i[H(t),p(t)} + £ %1Q[L k ]p(t), (2) 

k=\ 



where 



D[L k ]p = L k ph\ - l -L\L k p - l -pL\L k . 



For such a single qubit system, we can divide H(t) into three parts H(t) = Hq + H^ + H u , where 
the free Hamiltonian is Ho = j<j z , the control Hamiltonian is H u = Y i j= x , y ,z u j( t )^j' ( u j(t) £ R> 
7) = jCF,), and the uncertainties in the system Hamiltonian are = co(t)I z + £ x (t)I x + £ y (t)I y 
(G)(t),£ x (t),£y(t) e R). The Pauli matrices o = (o x ,o y ,o z ) take the following form: 

I ° 1 I 





i 

is the first class of uncertainties we will consider in this paper. The unitary errors in lETI 
belong to this class of uncertainties, and one-qubit gate errors also correspond to this class of 
uncertainties ll28l . A second class of uncertainties are uncertainties Sy k residing in the coupling 
strength y k . Since the Lindblad equation is an approximate equation for the open quantum system 
coupling with its environment, this class of uncertainties may come from inaccurate modeling 
as well as time-varying coupling between the system and environment. We assume that all the 



uncertainties are bounded, i.e., \(0{t)\ < co, j£^{t) + £^ < £ and \5y k \ < y, where constants 
CO > 0, £ > and 7 > are given. 

For a qubit system, its state p can be represented in terms of the Bloch vector r = (x,y,z) = 
(tx{po x },tr{po y },tr{po z }): 

I + ra 

P = —■ (4) 

After representing the state p with the Bloch vector, the pure states (i.e., with tr(p 2 ) = 1) for the 
qubit system lie on the surface of the Bloch sphere and the mixed states (i.e., with tr(p 2 ) < 1) 
occupy the interior of the Bloch sphere. The purity of p is defined as P = tr(p 2 ). A pure state 
can also be represented by a unit vector in a complex Hilbert space, where p = |y/}(i//|, 
= (\y)y and the operation refers to the adjoint of X. The fidelity of an arbitrary state 
p in terms of can be defined as (yf\p\y/). Thus, the fidelity between two pure states 
and |0) reduces to |(i//|0)| 2 . A projective measurement with o z on the qubit in state p will 
make the state collapse into |0) with probability (0|p|0) or into |1) with probability (l|p|l) 
(such a process is referred as the measurement collapse postulate), where |0) and |1) are the 
eigenstates of o? with corresponding eigenvalues 1 and -1, respectively. Another useful quantity 
is the coherence which can be defined as C = x 1 + y 2 , where x = tr(p<7 x ) and y = tr(po y ) (see, 



e.g., (581, Il59l , lfT3l ). A decoherence process due to the interaction of a quantum system with 
its environment may reduce its purity or coherence. 

We will consider the following four cases in this paper: 

A) No decoherence (i.e., Yk = 0). This case corresponds to a closed quantum system with a 
pure state p t satisfying the Schrodinger equation 

pt = -i[H(t),p t ]. (5) 

B) Amplitude damping decoherence. In this case, the population of the quantum system can 
change (e.g., through loss of energy by spontaneous emission). The evolution of p t can be 
described by the following equation: 

p t = -i[H(t),p t ]+y t (o-p t o + - ^o + 0-p t - ^p t o + oJ) (6) 

where o = j(<j x — io y ), c+ = ^(<J x + i(J y ) y t = Yo + 8y t and \8y t \ < J. We also assume that 
To > 7, which guarantees the coupling strength y > 0. 

C) Phase damping decoherence. In this case, a loss of quantum coherence can occur without 
loss of energy in the quantum system. The evolution of the state may be described by the 
following equation: 

p t = -i[H(t),p t ]+y t (o z p t o z -p t ). (7) 

D) Depolarizing decoherence. This decoherence maps pure states into mixed states. The 
dynamics can be described by the following equation: 

pt = -i[H(t),p t ] + Yt(o x p t o x - pt) + YtiOyptOy - p t ) + Yt(o z ptO z - p t ). (8) 

The objective of this paper is to design control laws for single qubit systems guaranteeing 
required robustness with the two classes of uncertainties. The required robustness for the four 
cases above is defined using the concept of a sliding mode domain, respectively, as follows. 

Definition 1: ll29ll The sliding mode domain for a single qubit system without decoherence 
(closed systems) is defined as % = {\w} '■ \ (^\¥)\ 2 > 1 -Po,0 < po < 1}. 

Definition 2: The sliding mode domain for an open qubit system with amplitude damping 
decoherence is defined as 2# a = {p ■ (0[p ]0) > 1 — po,0 < po < 1}. 

Definition 3: The sliding mode domain for an open qubit system with phase damping 
decoherence is defined as Si p = {p : x 2 +y 2 > C,x = tr(pa x ),y = tr(pa v ),0 < C < 1}. 



Definition 4: The sliding mode domain for an open qubit system with depolarizing deco- 
herence is defined as ^ = {p : trp 2 > P,0.5 < P < 1}. 

Remark 1: The definition of 9> c implies that the system's state has a probability of at most 
po (which we call the probability of failure) to collapse out of <2i c when making a projective 
measurement with the operator o z . We aim to drive and then maintain a single qubit's state in the 
sliding mode domain S) c . However, the uncertainties may take the system's state away from 
& c . The sampling process (a measurement operation) unavoidably makes the sampled system's 
state change. Thus, we expect that the control law will guarantee that the system's state remains 
in @ c , except that the sampling process may take it away from @ c with a small probability (not 
greater than po). The definition of S> a has a similar meaning to S> c . The difference lies in the fact 
that the quantum state in Definition [2] could be a mixed state p and the system is also subject to 
amplitude damping decoherence. From Definition |3l we know all states in Qi p have coherence 
of at least C. Definition |4] defines ^ as a set where the purity of an arbitrary quantum state is 
not less than P. 

III. Main Methods and Results 

In this paper, we propose a sampled-data design method for robust control of quantum systems 
with uncertainties. A key task is to design a sampling period as large as possible to guarantee 
the required robustness defined using a sliding mode domain. The sampling process is taken 
as an important control tool to modify the system dynamics nonunitarily. For the cases without 
decoherence and with amplitude damping decoherence, it is also necessary to design a control 
law to drive the system's state back to the corresponding sliding mode domain when the sampling 
process makes the system's state collapse out of the sliding mode domain. Such a control law 
corresponds to a unitary transformation and we refer to it as "unitary control" in this paper. 
The sequel will provide the main methods and results for the four cases of uncertain quantum 
systems and then present some illustrative examples. 

A. No decoherence 

The objective is to develop a control strategy to guarantee the required robustness when 
bounded uncertainties exist in the system Hamiltonian. According to Definition 1, we specify 
the required robustness as follows: (a) maintain the system's state in the sliding mode domain 



Ql c in which the system's state has a high fidelity (> 1 — po) with the sliding mode state |0), and 
(b) once the system's state collapses out of @ c upon making a measurement (sampling), drive it 
back to Qi c within a short time period /3T C and maintain the state in <2i c for the following time 
period (1 — /3)T C (where < /3 < 1 and T c is the sampling period). /3 is used to characterize the 
proportion of time that the unitary control is applied within the corresponding sampling period. 
Generally we choose /3 to satisfy <C 1, and this assumption will be helpful for designing the 
unitary control, which is demonstrated in the examples. To guarantee the required robustness, we 
design a control law based on sampled-data measurements as follows: For any sampling time nT c 
(n = 0, 1,2, . . .), (i) if the measurement data corresponds to |0), let the system evolve with zero 
control and sample again at the time (n+ l)T c ; (ii) otherwise, apply a unitary control to drive 
the system's state back into a subset S c of 3> c from the time nT c to (n + j3)T c , then sample again 
at time (n+ l)T c . The control operation is switched between (i) and (ii) based on the sampled 
data measurements. In (ii), to guarantee the desired goal when t £ [(n+ 1 — j5)T c , (n + 1)T C ], the 
unitary control should drive the system's state into S c C & c . S c can be defined as S c = {\\ff} : 
|(0|i//)| 2 > 1 — apo,0 < po< 1,0 < a < 1}. The sampling period T c and the unitary control can 
be designed offline in advance. The basic method we use is illustrated in Fig. CD The sequel 
will outline the design of the sampling period and establish a relationship between a and /3 to 
guarantee the required robustness. 

Using a similar argument to Theorem 1 in [29], we have the following result [|30l . 

Lemma 5: For a single qubit with initial state |y(0)) = |0) at time t = 0, the system evolves 
to under the action of H(t) = [1 + co(t)]I z + £ x (t)I x + £ y (t)I y (where |co(/)| < ft), ft) > 0, 



the state will remain in % = {\yr) : \ (0\\ff)\ 2 > 1 — po} (where < po < 1). When a projective 
measurement is made with the operator ov at time t, the probability of failure p = |(1| !//(?)) | 2 
is not greater than pq. 

We use T c defined in © as the sampling period to guarantee the required performance. If the 
sampling data corresponds to |1), a unitary control is required to drive the state back to a subset 
$ c of 2> c . The following theorem gives a sufficient condition on the relationships between a, po 
and /3 to guarantee the required robustness. 




£%(t) +£, 2 (0 < £ and e > 0). If t e [0,T C ], where 



arccos(l —2po) 



(9) 



£ 




Fig. 1. The proposed sampled-data control scheme for a single qubit system without decoherence. The labels "Control", 
"Sampling" and "Uncertainties" refer to the evolution process under the unitary control, the sampling process, and uncertainties 
in the system Hamiltonian, respectively. 



Theorem 6: For a single qubit with initial state |i/f(0)) satisfying |(i//(0)|l)| 2 < apo (0 < 
a < 1) at time t = 0, the system evolves to \y/(t)) under the action of H(t) = [1 + G)(t)]I z + 
e x (t)I x + e y (t)Iy (where yj e%{t) + e%{t) < e, e > 0, |a>(f)| < co and co > 0). If t e [0, (1 -j8)r c ] 
and 

a< l-cos[/3arccos(l-2p )] 
2p 



where < /3 < 1 and 



_ arccos(l - 2gg) 
7c " e ' (U) 



then the state will remain in ^ c = {\\j/) : | (0|y) | 2 > 1 — po} (where < po < 1). When a projective 
measurement is made with the operator a, at time t, the probability of failure p = |(1| v(0) | 2 
is not greater than pq. 

Remark 2: Using Lemma |5] and Theorem |6] we aim to maintain the state in $ c by imple- 
menting periodic sampling with period T c in ©. This theorem provides a sufficient condition to 
guarantee the required robustness. Given po, /3, we can select a satisfying (flOl) in Theorem |6] 
If the sampled result is |1), we apply a unitary control to drive the state into S c . The sampling 
period and the unitary control can be designed in advance. Different approaches can be used to 



design such a unitary control law. In this paper, we will employ a Lyapunov method Il60ll - ll63l 
in Example 2 to accomplish this task for the closed quantum system. 

Remark 3: The design scheme above involves a sampling process and a unitary control. It 
is similar to the approach used in [|2"9"1 . The difference lies in the fact that the scheme in this 
paper involves a fixed sampling period T c . However, the approach in [1291 involves at least two 
measurement periods T (equivalent to T c in this paper) and T\ {T\ T). This situation means 
that the approach of [1291 may require measurements which are very close together, which may 
be difficult to achieve in practice. In this sense, the sampled-data design in this paper is more 
practical than the method in [|2~9~1 . 



B. Amplitude damping decoherence 

For single qubit systems with amplitude damping decoherence, if the initial state is excited 
state |0), the decoherence will drive this excited state to the ground state |1). The objective is 
to design a control law to guarantee the required robustness defined by @ a . We use a similar 
sampled-data design method to that in the case without decoherence. That is, if the state is |0) 
at t = nT a (n = 0, 1,2, . . .), we design a sampling period to maintain the system's state in *3) a 
by implementing periodic sampling with period T a ; if the measurement makes the state collapse 
into |1) (with a probability p < po), we design a unitary control to drive the state back into a 
subset S a of $> a from t = nT a to (n + j3)T a , and then sample again at t = (n+ l)T a . In order to 
determine the required sampling period, we have the following results. 

Theorem 7: For a single qubit with initial state |0) at time t = 0, the system evolves to p(t) 



subject to © where H(t) = [1 + a>(t)]I z + e x (t)I x + e y {t)I y (Je%(t) + e*(t) < e, £ > 0, \(o(t) \ < 
CO and to > 0) and the coupling strength of amplitude damping decoherence is y t = yo + 8y t 
(\8Yt\<Y)- Kte[0,T a ] with 

T a = 2P0 , (12) 

v/4e 2 + (7o + r) 2 + (7o + r) 

the state will remain in £t a = {p : (0|p|0) > 1 — po} (where < po < 1). When a projective 

measurement is made with the operator o z at time t, the probability of failure p= (l|p|l) is not 

greater than pq. 

Corollary 8: If po < j — — j= ^ y = ^ 2 and t G [0,7^], the sampling period T a in (TT2l) can 



be replaced by T' a to guarantee the same robustness as in Theorem |7J where 

T a = 2po . (13) 

4eJp -pl + 2(ro + r)^-Po) 



When there exist no uncertainties in the system Hamiltonian (i.e., = 0), the sampling period 
can be designed using the following proposition. 

Proposition 9: For a single qubit with initial state |0) at time t = 0, the system evolves to 
p(t) subject to © where H(t) = L and the coupling strength of amplitude damping decoherence 
is 7 t = 7o + 8y t (\8y t \ < y). If t e [0, T%] with 

ln(l- P0 ) 

To + 7 ' ( ] 

the state will remain in *2i a = {p : (0|p|0) > 1 — po} (where < po < 1). When a projective 

measurement is made with the operator o z at time t, the probability of failure p = (l|p|l) is not 

greater than pq. 

Remark 4: From the proof of Proposition it is clear that = — ln ^ {)+ ^ exactly corre- 
sponds to the case 8y t = y when = 0. In this sense, the sampling period T% = ~ 1 "^ ) +y°" > is 
optimal to guarantee the required robustness. From the proofs of Theorem [7] and Corollary |8] it 
is clear that T' a > T a . The relationship T" > T a for arbitrary po can be proved by the following 
steps: (a) Define F(p ) = T" - T a ; (b) observe F(p = 0) = 0; and (c) verify ^ > 0. 

Hence, for different situations we may use T a , T' a or T" as the sampling period to guarantee 
the required performance. If the sampled data corresponds to |1), a unitary control is required 
to drive the state back to a subset <S a of 3> a - The subset S a may be defined as S a = {p '■ 
(0|p|0) > 1 — ocpo,0 < po < 1,0 < a < 1}. The following theorem gives a sufficient condition 
on the relationships between a, po and /3 to guarantee the required robustness (The following 
conclusion is also true when T a can be replaced by T' a or T"). 

Theorem 10: For a single qubit with initial state p(0) satisfying (l|p(0)|l) < ocpo (0 < 
a < 1) at time t = 0, the system evolves to p(t) subject to © where H{t) = [1 + Co(t)]I z + 



£ x (t)I x + £y(t)I y (J£l(t) + ej(t) < e, e > 0, \co(t)\ < co and co > 0) and the coupling strength of 
amplitude damping decoherence is y t = 70 + ^7 (\8y t \ < 7). If t G [0, (1 — /3)rj and 

a</3 (15) 

where < /3 < 1 and 

T - 2po 



the state will remain in 3 a = {p : (0|p|0) > 1 — po} (where < po < 1). When a projective 
measurement is made with the operator ov at time t, the probability of failure p= (l|p|l) is not 
greater than pq. 

C. Phase damping decoherence 

For a single qubit with phase damping decoherence, we define the coherence as C = x 2 +y 2 
where x = tr(pc x ) and y = tr(p<7 y ). The phase damping decoherence will reduce the coherence 
of the system. The objective is to guarantee that the state has coherence not less than C by 
periodic sampling when there exist uncertainties in the coupling strength of system-environment 
interaction and in the system Hamiltonian. To determine the required sampling period, we have 
the following results. 

Theorem 11: For a single qubit with initial state po satisfying Co = [tr(poC7 Y )] 2 + [tr(po<7 y )] 2 = 
1 at time t = 0, the system evolves to p t subject to © where H(t) = [1 + co(t)]I z + e x (t)I x + £ y (t)I y 



(Je 2 (t) + £y(t) < £, £ > 0, \co(t)\ < co and co > 0) and the coupling strength of the phase 
damping decoherence is y = Yo + dy t (\8y t \ < y). If t e [0,T p ] with 



1 — C „,V,q„ Mrs, I r,,\2 \ c.2. 



T = { 4 ^2(ro+r)' (17) 
" '-^P^, when 4(70 + 7)2 < e\ 

the state will remain in @ p = {p t : [tr(p t o x )} 2 + [tr(p t o y )] 2 > C,0 < C < 1}. When a periodic 
projective measurement is made with the operator o x on the system, the sampling (measurement) 
period T p can guarantee that the state remains in @ p . 

Corollary 12: When £ 2 = 2(7 + 7) 2 , the sampling period T p in (IT7T) can be replaced by T' p 
to guarantee the same robustness as in Theorem [TT] where 

T t = 1 ~ ^£ (18) 

p 2^(70 + 7)' 

If = 0, we can design the sampling period using the following proposition. 

Proposition 13: For a single qubit with initial state po satisfying Cq = [tr(po<7 x )] 2 + [tr(po<7y)] 2 
1 at time t = 0, the system evolves to p (t) subject to © where H(t) = I z and the coupling strength 
of the phase damping decoherence is y = Jo + 8y t (|<5y| < 7). If t e [0,T' r ] with 



when 4(70 + y) > e > 



the state will remain in @ p = {p t : [tr(p t a x )] 2 + [tr(p t o y )] 2 > C, < C < 1 }. If a periodic projective 
measurement is made with the operator o x , the sampling period Tp can guarantee that the 
system's state remains in @ p . 

Remark 5: For 2(y+yo) 2 = £ 2 , it is straightforward to prove that T' p > T p . The relationship 
Tp > T p can be proved by the following steps: (a) Denote Y = 1 — C; (b) define F(Y) = Tp — T p ; 
(c) observe F(Y = 0) = 0; and (d) verify > 0. From the proof of Proposition [T3l it is clear 
that the sampling period Tp is optimal to guarantee the required robustness when = 0. For 
this case with phase damping decoherence, we can also make projective measurements with the 
operator a y , which does not affect the conclusions. Moreover, in this case, no unitary control 
is required and measurement is the only tool needed for guaranteeing the required robustness. 
It is worth mentioning that several methods based only on measurements have recently been 
proposed for controlling quantum systems (see, e.g., [1641 - 116710 . 

Remark 6: We can also consider a class of imperfect measurements. This class of uncer- 
tainties may arise from precision limitations of the measurement apparatus or from system errors 
in the measurement device. Measurement with the operator o z will make the system collapse 
into |0) or 1 1 ) (eigenstates of Hq). We consider the imperfect measurement model shown as in 
Fig. [2l poi is the error probability of measurement from |0) to |1), that is, the probability that 
one obtains the result |1) when making a measurement on the system in |0); pio is the error 
probability of measurement from |1) to |0), where < /?oi < 1 an d < pio < 1. This class of 
imperfect measurements does not affect the effectiveness of the sampled-data design. Thus, the 
proposed method can tolerate this additional uncertainty in the sampling process. 

D. Depolarizing decoherence 

For a single qubit, the depolarizing decoherence will reduce the purity P = tr(p 2 ) of the 
system's state. The objective is to guarantee that the purity of the state is not less than P by 
periodic sampling when there exist uncertainties in the coupling strength of system-environment 
interaction and in the system Hamiltonian. To determine the required sampling period, we have 
the following results. 

Theorem 14: For a single qubit with initial state po satisfying tr(pp) = 1 at time t = 0, the 
system evolves to p t subject to ® where H(t) = [l + co(t)}I z + £ x (t)I x + £y(t)I y {^J £ 2 {t) + £ 2 (t) < 
£, £ > 0, \co(t) \ < (O and 0) > 0) and the coupling strength of depolarizing decoherence is 
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Fig. 2. The model for imperfect measurement, poi is the error probability of measurement from |0) to |1) and pio is the error 
probability of measurement from |1) to |0). 



Yt = Yo + 8r, (\Sy t \ < Y). If t e [0,T d ] with 

ln(2P-l) 

r '=-»W' (20) 

the state will remain in ^ = {p t : tr(p f 2 ) > P, 0.5 <P < 1}. If periodic projective measurements 
are made with the operator o z , the sampling period Td can guarantee that the state remains in 

Remark 7: The selection of measurement operators (i.e., o x , o y or a z ), uncertainties in the 
system Hamiltonian ^ or = 0) and the imperfect measurement described in Fig. |2]do 
not affect the conclusion in Theorem [14] The sampling period Td is also optimal to guarantee 
the required robustness. 

The sampling periods for the different cases considered above are summarized in Table 1. 



Table 1: Summary of sampling periods for different cases. H A (x,y,z) = co(t)I z + £ x (t)I x + 
£y(t)I y (where J 1 £%(t) + ej(t) < e, £>0, \co(t) \ < co and go > 0), /(e,y ,y) = \ 



2 v / 4e 2 +(ro+r) 2 ' 

and coupling strength y = yo + <5y (|<5y| < 7)- H A (x,y,z) is also considered for the two cases 
Po < /( £ 5 7o,7) and £ 2 = 2(y + y) 2 . The parameters values p = 0.01, Jq = 0.9, 7= 0.1 and 
C = P = 0.95 are assumed for the calculation of the right two columns. When 4(yo + y) 2 > e 2 , 



cases 



sampling period 



0.2 



V2 



closed system with H^(x,y,z) 



arccos(l— 2po) 



1.0017 



r c = 0.1417 



amplitude 
damping 
decoherence 



HA(x,y,z) 



T a 



2po 



T a = 0.0096 



T a = 0.0050 



Po<f(£,Yo,Y) 



r 



2po 



4e v /po-p ( 2 ) +2(yo+y) ( 1 - P0 ) 



T a = 0.0079 



H A = 



rj.ii _ ln(l-po) 

yo+y 



T' a ' = 0.0101 



T' a ' = 0.0101 



phase 
damping 
decoherence 



2(7o + rf 



H A = 



i-Vc 



T' - 



P 



InC 



4(yo+y) 



T p = 0.0088 



T p = 0.0088 



Tp = 0.0090 



Tp' - 0.0128 



T p ' = 0.0128 



depolarizing 
decoherence 



H A = 



ln(2P-l) 

' 3(%+r) 



ln(2P-l) 

' 8(ro+y) 



T d = 0.0131 



r rf = 0.0131 



0.0131 



r d = 0.0131 



E. Illustrative examples 

Example 1 (Sampling periods): The values of sampling periods are shown in the right two 
columns of Table 1 for several specific cases, where we have assumed po = 0.01, jo = 0-9, 
y = 0. 1 and C = P = 0.95. Further, we can consider a real quantum system of a superconducting 
box in 0, ll68l . Let the resonance frequency Q)q = 2% x 100MHz and the cavity decay rate % — 
2% x 0.8MHz. Assume that f— 2?rxQ 9 8MHz and £ = 2nx 1.0MHz. Hence, the cavity decay time 
Ty= 198.9ns. Using the results in Table 1, we can get the real sampling periods as T c = 8.0ns, 
f a = 1.0ns, t' a ' = 1.8ns, T p = 1.6ns, Tp = 2.3ns and Tj = 2.5ns. 

Example 2 (Unitary control for Case A)): Theorem |6] gives a sufficient condition for de- 
signing a unitary control to guarantee the required robustness. Here we employ a Lyapunov 
method Il60l - [l62l to design such a unitary control where the Lyapunov function is constructed 



based on the Hilbert- Schmidt distance between a state \\jf) and the sliding mode state |0); i.e., 
V(|y), |0)) = j(l - |(0|v/)| 2 ). The control values can be selected as (for details, see [|29), ll63l): 

H , = ^A(3[e'' z ^^(0|4|vA)]), (* = x,y,z) (21) 

where 3[a + bi] = b (a,b G R). Here Zc denotes the argument of a complex number c, the 
parameter K^> may be used to adjust the control amplitude, and /(•) satisfies xf(x) > 0. We 
define Z{i//|0) = 0° when (v|0) = and adopt the parameter values of po = 0.01, £ = 0.2 and 
/3 = 0.05. From the simulation in ||29l , we find that the Lyapunov control is not sensitive to small 
uncertainties in the system Hamiltonian. Additional simulation results suggest that the robustness 
of the Lyapunov control can be enhanced if we choose the terminal condition |(l|i//(?))| 2 <r\apQ 
(where < r\ < 1) instead of |(l|ty/(?))| 2 < apo. Here, we select 77 =0.8. Hence, we design the 
sampling period T c = 1.0017 using ©. Using Theorem |6l we select a = 2.5 x 10 -3 . We design 
the Lyapunov control using (l2TT ) and the terminal condition | (1 1 i//(f)}| 2 <rj(Xpo = 2.0 x 10~ 5 with 
the control Hamiltonian H u = \u{t)o y . Using we select u(t) = K(Z[e iA ^^ (0\o y \y{t))\), 
K = 500, and let the time stepsize be 8t = 10 . We obtain the probability curve for |0) shown 
in Fig. [3£a) and the control value shown in Fig. Hfb). For the noise £(t)I x or e(t)I y where e(t) 
is a uniform distribution in [—0.2,0.2], additional simulation results show that the state is also 
driven into S c using the Lyapunov control in Fig. |3{b). 

Example 3 (Unitary control for Case B)): For amplitude damping decoherence, Theorem 
[l0]gives a sufficient condition for designing a unitary control to guarantee the required robustness. 
Here, we employ a constant control H u = luo y (u = 6466), and assume po = 0.01, £ = 0.2, 
70 = 0.9 and /3 = 0.05. Using Theorem \\0[ we may select a = 0.05. Let the time stepsize be 
8t = 10 -7 . The curve for zt = tr(p t o z ) is shown in Fig. BJ 

Remark 8: The required unitary control can be designed using strategies such as the Lya- 
punov method and optimal control theory. In Example 2 and Example 3, we used simulation to 
find appropriate control amplitudes for achieving the required objectives. Additional simulation 
experiments also show that u(t) can tolerate small uncertainties. Since it is necessary to drive 
the state back to a subset of the sliding mode domain within a small time period, the required 
control amplitudes generally are relatively large, which is similar to the case of decoherence 
control based on dynamical decoupling [|24l. [PT71 . The selection of small /3 makes it reasonable 
that we first design the unitary control by ignoring possible uncertainties and then verify the 
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Fig. 3. The probability of |0) and the control value u(t) under the Lyapunov control. 
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Fig. 4. The curve of zt for the case with amplitude damping decoherence, where H u = 3233(7 V . 

robustness of the unitary control to uncertainties by simulation. Here, we present only simulated 
examples to demonstrate how such a unitary control can be designed. A systematic investigation 
into the design of the unitary control and finding optimal control amplitudes that can tolerate 
uncertainties will be the subject of future work. 



IV. Proof of the Main Results 

A. Proof of Theorem |6| 

To prove Theorem |6] we first prove two lemmas (Lemma [15] and Lemma [T6l) . Lemma 
compares the probabilities of failure for H{t) = [1 + (o{t)]I z + £cos0 o A> + £sin0 o / T and H{t) = 
ecos^oA' + £sm ^o4- Lemma [TBI together with Lemma [T6l demonstrates that H = el x can be used 
to estimate an upper bound on the probability of failure for H(t) = [1 + 0)(t)]I z + £ x (t)I x + £ y (t)I y 
when zo = 1 . 

Lemma 15: For a single qubit with initial state (xo,yo,Zo) = 1), the system evolves 
to (xf,yf,zf) and (xf,yf,zf) under the action of H A = [l + 0)(t)]l z + £cos^ol y + £sm(j>ol x (with 
constant e > and \co(t)\ < (0) and H B = e cos fol y + e sin 00^-, respectively. For arbitrary f e 



[0, 



2\/4+i 



Proof: For the system with Hamiltonian // A = [1 + 0)(t)]I z + £ cos ^oI y + £sin(j>oI x , using © 
and ©, we obtain the following state equations 



/ 



V 



-[1 + ©(/)] ecos^o 
l + co(/) -esin^o 

-£cos0n £sin0n 



\ 



( + \ 



4 

\4 J 

where (4,$, 4) = (0,0,1)- 

Consider (o{t) as a control input and select the performance measure as 



/ 



4 
4 



(22) 



(23) 



\ ( x \ 



J(00) = Zf. 

We introduce the Lagrange multiplier vector X(t) = fait), ^(t)) 7 and obtain the corre- 

sponding Hamiltonian function as follows: 

^ -[1 + co(t)} £COS0 O 
m(r(t),(0(t),l(t),t)=l T (t) l + a>(t) -£sin0 o 

^ — £cos0o £ sin 0o 
where r(t) = (x t ,y t ,Zt)- That is 
U(r(t),(o(t),X(t),t) 

= [1 + ft)(0](A 2 (0^ - + £cos0o(Ai(Oz f - h(t)x t ) - esin0o(A 2 (?)^ - ^(O^)- 

(25) 



/ 



37 



(24) 



According to Pontryagin's minimum principle [1691 , a necessary condition for (0*(t) to minimize 
/(ft)) is 

H(r* (t), co* (t) , V (t) , t) < H(r* (t) , <o(f) , V{t) ,t) . (26) 

Hence, if we do not consider singular cases (i.e., Xi{t)x t — X\{t)y t = 0), the optimal control o)*(f) 
should be chosen as follows: 



(0*{t) = -cosga(X 2 (t)xt-Xi(t)yt). 



(27) 



That is, the optimal control strategy for G)(t) is bang-bang control; i.e., (Q*(t) = ft) = +0) or — ft). 
Now we consider H A = (1 + cb)I z + e cosily + £sin0o/ x , which leads to the following state 
equations 

' \ ( + £cos0 o \ ( x$ ^ 

l + ft> — £ sin 0o 

-e cos 0o £ sin 0o 



yf 



\ 



) 



4 



(28) 



where (*o,)>g,Zo) = (0,0, 1). The corresponding solution is 

/ 



(4) 




yf 




U J 





ecos ( 



sin w _ (^fft cos W + 

a/(1+<b) 2 +£ 2 (1+co) 2 +£ 2 (l+ffl) 2 +£ 2 

esindm (l+a)ecos0o , (l+d>)ecos( 

(l+6) 2 +£ 2 



sin W 



: cosW + — 



(l + ffl) 2 +£ 2 



V 



(l+<3) 2 +£ 2 C0S Vt + (l+i"+e 2 



(29) 



/ 



where v = a/(1 + ft)) 2 + £ 2 . From (|29l , we know that z t is a monotonically decreasing function 
in t when t G [0, 2 ^* +£2 ] • Hence, we only consider the case ? G [0,?/] where *y G [0, — /J=f ]■ 

Now consider the optimal control problem with a fixed final time ?y and a free final state 
Yf = (xf,y/,Zf). According to Pontryagin's minimum principle, X*(tf) = ^r*(tf), and it is 
straightforward to verify that ^(tf), X^(tf)) = (0,0, 1). Now consider another necessary 



condition X (t) 



k{t) 



dr 




' k{t) N 













which leads to the following relationships: 





1 + ft) 

-£cos <t 



+ 


£ sin 0o 



£ cos 0o 
— £ sin 0o 




^ ( m ) 



(30) 



where {X\{tf), ^(f/), X$(tf)) = (0,0,1). The corresponding solution is 



f MO 












{ m 


) 





( 



tuu " TO -in riff- t\ v 1 +^) £sin( fe rcKXiiu t) I ( — 



esin( 



= sinu(f/- 



[l+©) 2 +£ 2 



,l+<B) 2 +e 2 



[+ffl) 2 +e 2 



1+©) 2 



(3D 



We obtain 



X 2 {t)x t - X\{t)y t 



e 2 (l + <y) 

3/2 



V 



[sin vt + sin v(tf — t) — smvtf]. 



(32) 



It is easy to show that the quantity (X2(t)x t — X\(t)y t ) > occurring in (127b does not change 



sign when tf G [0, 



2V4+i 



and f G [0,f/]. Hence, the optimal control is 5*(f) = (Q = -co. 



We now exclude the possibility that there exists a singular case. Suppose that there exists a 
singular interval [to,h] (where fo > and we assume that [fo,fi] is the first singular interval) such 
that when f G [to,t\] 

h{t) = X 2 {t)x t -X l {t)y t = Q. 



We also have the following relationship 

h(t) = X 3 (t)x t -X l (t) Zt = 
where we have used (1221 ) and the following costate equation 

( ki(t) \ I o 



(33) 



(34) 



k{t) 
\ hit) J 



\ 



J 



h{t) 
X 2 {t) 



(35) 



-[1 + ©(/)] ecos^o 
X(t) = X^(t) 1 + 0(0 -£sin0 o 

— e cos 0o e sin 0o 
If to = 0, we have (xo,yo,Zo) = (0,0, 1). By the principle of optimality (691, we may consider 
the case tf = t\. Using (|33T> , (|341 and ( Ai (?i ) , A2 (?i ) , A3 (?i ) ) = (0,0,1), we have x h = and 
y tl = 0. Using the relationship of xj + y 2 + zj = 1, we obtain z ?1 = 1 or — 1. If z fl = 1, the initial 
and final states are the same state |0). However, if we use the control 00(f) = ft), from (l29l we 
have Zt x {(b) = ( 1+ J)2 +£ 2 cosufi + ^^y.\_ £ i < Zt x = l- Hence, this contradicts the fact that we 
are considering the optimal case minzy-. If Zt x = — 1, there exists < t\ < t\ such that zi x = 0. 
By the principle of optimality 11691 . we may consider the case tf = t\. From the two equations 
(1331) and (|34|) . we know that zj = 1 which contradicts z~ tl = 0. Hence, no singular condition can 
exist if fg = 0. 



If to > 0, using (|271) we must select (o(t) = ft) when t G [0,/o]. From (|32l , we know that there 
exist no to G (0,tf) satisfying hi{tQ)x t0 ~~ hi(to)y tQ = 0- Hence, there exist no singular cases for 
our problem. From the previous analysis, co(t) = -co is the optimal control when t G [0, — 7=^]. 

For the system with Hamiltonian H B = £cos0o/y + 8sin0o/x> using © and ©, we obtain the 
following state equations 



yf 



( 



\ 



e cos (j>o 

— £ sin 0o 

-£cos0o £ sin 0o 



yf 



(36) 



where (x^,Vq,Zq) = (0,0, 1). The corresponding solution is 



yf 

\zf ) 

We define F(t) and f(t) as follows: 



( 



(l-£0) 2 + £ 2 



cos 0o sin et 
-sin 0o sin et 
coset 



(1-fti) 2 
cos vt + — — — cos et, 



(37) 



m=F(t) = - 

Now, consider t G [0, — 7=f ] to obtain 

/(0 



v /(l_ w) 2 + £ 2 



(l-£0) 2 + £ 2 

sinw + £sin£f. 



£ (cosef — cosuf) > 0. 



(38) 
(39) 

(40) 



It is clear that f(t) = only when t = 0. Hence f(t) is a monotonically increasing function and 



min/(0=/(0) = 0. 



Hence, we have 



(41) 



/(0 > 0. 

From this result, it is clear that F(t) is a monotonically increasing function and 

minF(f) =F(0) =0. 

Hence F(t) > when ? G [0, ^~/|=f ]• Therefore, we can conclude that > zf for arbitrary 
We now present another lemma. 



Lemma 16: For a single qubit with initial state (jco, yo, Zo) = (0, 0, 1), suppose that the system 
evolves to (x t ,y t ,Zt) under the action of H = £ (cos <j>I y + sin <j>I x ) (0 is a constant). Then, z? is 
independent of 0. 

Proof: For // = G)(sin0/ x + cos0/y), from (1371 ), we have 

z f = coser. 

It is clear that zt is independent of 0. ■ 
Remark 9: Since z* is independent of 0, it is enough to consider a special case = f when 

analyzing zr under H = £(cos(j)Iy + sin (j)I x ). 
Now we can prove Theorem |6] 

Proof: For a single qubit, assume that the state at time t is p t . If we make a measurement 

with the operator o z , the probability 77 that the state will collapse into |1) (the probability of 

failure) is 



p = (Mpt\i) 



l-z t 



For a closed single qubit system, its state | Xjf) can be represented as 

|V/)=cos||0) + ^sin||l), 



(42) 



(43) 



where its bloch vector corresponds to (x,y,z) = (sin cos <p, sin sin <p, cos 0), 6 [0, ft], <p G 
[0,2rc]. 

FOr // A = [1 + Q)(t)]l z + £ x (t)l x +£y(t)ly, using p = — i[H A ,p] and ©, we obtain the following 

state equations 



I -[l + fl»(0] £y(0 \ 



1 + Q)(f) 

^ -e y (t) e x (t) 







7 



(44) 



Define e(t) = Je}(t) + e*(t) and e x (t) = e(?)sin0, , £,,.(?) = e(f)cos0 f . This leads to the fol- 
lowing equation 

( &\ I -\l + 6)(t)} e(t)cosdk \ I 



yf 



-[1 + ©(/)] £(f)cos0 r 

1 + Q)(0 -e(?)sin0f 

-e(f)cos0f e(f)sin0 ; 



(45) 



( 



\ 



\ 



J 



yf 



(46) 



For H B = el x , we have 


-£ 
e 

When (jcq,3^,Zo) = O^o^o^o) = f° r ^ ~~ > ^> we nave ^ rom Lemma [T51 and Lemma [TBI 

4>4- (4V) 

We will now prove that the relationship zff > zf f (At — > 0) is also true for zf — Zq = cos 0o 
(where 0q G (0, tt)). We assume that there exist f G (0,Af] such that 



(48) 



Define f(t) = zf — zf- Since f(t) is continuous in t and /(0) = 0, there exists a time f* = 
sup{f|0 < t < t,f(t) = 0} satisfying /(**) = and /(/) < for t G (?*,?]. Hence 

f(t)\t=t*<0. (49) 

Let z$ = zf* =cosG*. We can assume (xf*,yf*,zf*) = (sin0*cos<p^sin0*sin<p*,cos0*) and 
(xf«,y?*,zf) = (O,-sin0*,cos0*) (where (p* G [0,2w]). Define A^(f) = -e(f) sin0*cos(0f + <p*). 
From (1451) and (l46l ). we have 



/(*)!*=** 

For arbitrary it is clear that 



zf\t=t*-zf\ t =t* = HmN(t) + £sm9*. 

t-tt* 



(50) 



N(t) > -esine*. (51) 

When iV(?) > -£sin0*, f(t)\ t = t * > 0, which contradicts (|49l . When iV(/) = -£sin0*, since 
sin0* 7^ 0, we have e(/*) = £ and ty t * = 27T — 9*. Using Pontryagin's minimum principle [|69l 
and a similar argument in Lemma [15] and Lemma dU we can prove zf* , ^ > z^ + at Hence we 
can conclude that for z$ = Zq = cos ^0 (where 0q G [0, it)) and Af — > 0, 

(52) 

From (BBl , we know that zf = cos(0 o + £O- When < t < ^r 1 , zf decreases monotonically 
in t. We now define g(t) = zf — zf and assume that there exist t = t\ G [0, 7r ~ e ° ) such that 
zf < zf r That is, g(t\) < 0. Since is continuous in t and g(0) = 0, there exists a time 

= sup{f|0 < t < t\,g(t) = 0} satisfying g(t) < for t G However, we have established 



that for any zf = zf and At ->■ 0, zf +At > zf +Ar which contradicts g(t) < for t e (t*,t\]. Hence, 
we have the following relationship for t G [0, 

4>i- (53) 

From (1421 ). it is clear that the probabilities of failure satisfy pf = — ^ < pf = ~^f~- That is, the 
probability of failure pf is not greater than pf for t E [0, ■^^ a ). 

Since zf = cos(6o + £t), we have Az| r = cos 0o — cos(0o + ej8r), where 

= arccos(l-2p ) 

When |(ty/(0)|l)| 2 < ocpo> using the previous argument, we have 

zi Tc > l-2apo + cos(0o + ef}T c ) -cos0 o = Af. 

Now let 

P=^— < — </>(). 

Using the fact 6o = arccos(l — 2apo), we have the following relationship 

a < l-cos[(l-/3)arccos(l-2p )] 
2p 



5. Proof of Theorem 

Proof: For the open qubit system subject to ©, when H(t) = [I + 0)(t)]I z + e x (t)I x + e y (t)I y 
(Je^(t) + e^(t) < £, £ > 0, \(o(t) \ < (0 and co > 0), y = yo + <5y (\8y t \ < y), using ©, we have 

\ 



{ * \ ( -Uro + 8y t ) -(l + «)(0) 



ft 



1 + qj(0 -j(7b + 5y t ) 



yt 



£y(t) 

xu.. , ^ _ £x(t) 
y -£,.(0 e x (0 -(y + 5y) y 

where (xo,;yo,Zo) = (0,0, 1). From (|56l , we have 

z r = -%(?)^ + e x (t)y t - (y + 8y t ) [z t + 1) 
>-2eVl-z?-(% + y)(zr + l)- 

Denoting 



+ 







^ -(yo + Sy) y 



(56) 



/(z)=2£Jl-z f 2 + (y + y)(l+z f ), 



(57) 



(58) 



we have 



Let = to find the solution z = . J 0+r Hence, 



d/(z) (ro+r)-2 £ ^. (59) 



max/(z) = /( Ia J + Y =) = J4e* + ( 7o + 7) 2 + (7o + 7)- (60) 



V4e 2 + (7o + 7) 2 

Hence, 



z f > -max/(z) = -y / 4e 2 + (ro + r) 2 -(ro + 7)- (61) 

When f G [0, r a ] where 

T a = == , (62) 

v/4£ 2 + (7o + 7) 2 + (7o + r) 

we have 

zr > 1 - (max/(z)> > l-2p . (63) 

Therefore, if one makes a measurement on the system with a z , the probability of failure (1 \p t | 1) = 
^<Po. ■ 



Proof: When po < \ — — ^ ^ y =^, from the proof of Theorem [7] we know for z G 



C. Proof of Corollary |#| 

Proof: 
[l-2po,l], 

max/(z) = f(l - 2p ) = 4e ^/p -p§ + 2(% + 7) ( 1 - Po) ■ (64) 
Hence, if t G [0, Tj] where 

rpl 2^0 

T a = 1 j (65) 



4e * //>o - + 2(7o + 7) ( 1 - A)) 



z t > I - [4e^/ po- pl + 2(y + y)(l - p )]t > I -2p . (66) 
It is clear that the probability of failure (l|p f |l) < pq. ■ 



D. Proof of Proposition [9] 

Proof: When H(t) = I z and y t = 70 + 8y r , the state equation of the system in © is 



( 



4(7o + $Yt) 



1 







\ 



1 4(7o + 8y t ) 



V 



7 



yt 



( 



+ 



+ 

where (xo,yo,zo) = (0,0, 1). It is clear that 

z t = -(Yb + 5y t )(l+z t ) > -(7o + 7)(l+^) 

From (1681) . we have 

Zt >2e-^+r)t_ h 
_„ ln(l-po) 








If * e [0, r a "] where 



7o + 7 



we have 



That is, the probability of failure (1 \p t | 1) < pq. 



E. Proof of Theorem \ 

Proof: From the proof of Theorem [7] we know 



z t > -max/(z) = -y / 4e 2 + (7o + 7) 2 -(7o + r). 
Now if the initial state zo > 1 — 2ocpo and ? e [0, (1 — J3)r a ], the system's state satisfies 



z f >zo-(V4£ 2 + (7 + 7) 2 + (7 + 7))(l-/3)r fl > 1-2(1 + a-p> . 
When a < fi, we have the following relationship 

Zf > 1 - 2p - 



Hence, the probability of failure satisfies (l|p f |l) < po. 



F. Proof of Theorem [771 

Proof: For a single qubit subject to © when H(t) = [\ + CO(t)]I z + £ x (t)I x + £ y (t)I y (\co(t)\ < 
co, Je^(t)+e^(t) < e, co > and £ > 0), y t = y + 8y t (\8y t \ < y), using ©, we have 



it 



( -2(y + 8y t ) -(l + ffl(0) %W \ / 



V 



l + co(t) -2(y + 8y t ) -£ x (t) 

-£y(t) £ x (t) J 

where Co = x\ + Vq = 1. Let C = xf + y 2 . We have 

C, = 2x t x t + 2y t y t = -4(y + 5%) (x 2 +y f 2 ) + 2z t (£y(t)x t - £ x {t)y t ) 
>-4(y + 8y t )(xj+yj)^2£y/l-(xi+yf)(\x t \ + \y t \) 
> _ 4 ( 7o + 5y,)(^ + y j) - 2£^\-{xf+yf)^2{xf+ y f) 

Let JV f = 2(7q + y) 2 C 2 - £ 2 C 2 + e 2 Q. We have 



(74) 



(75) 



2(y + 5y)C + £ V / 2C(l-C) < V2~V 4 (tt + 5 7/) 2 C? + 2e 2 Q(l -Q) < 2^ (76) 



Hence, Q > — A^JmaxN t . According to the definition of N t , it is easy to verify the fact 

2(7o + y) 2 , when 4(y + y) 2 > e 2 ; 



max N t 



(77) 



4e^-8(ro+7) 2 



when 4(y + y) 2 < £ 2 



If * G [0, r p ] where 



,2. 



jW^VV when4(y + y) 2 >^ 
P | (i z c)Vp^±g ; wh en4(y + y) 2 < £ 2 , 

we have C > C. 

G. Proof of Corollary U~2\ 

Proof: When e 2 = 2(y + y) 2 , from (EU and d76l we have 

C > -4 v / e 2 Q. 
It is easy to obtain the following relationship 



(78) 



(79) 



2dVC r > -4edf , 



(80) 



VQ > l -let. 



(81) 



If * G [0, T'] where 



we have Q > C. 



p 2^(70 +r) : 



if. Proof of Proposition 1731 

Proo/: When i/(f ) = 7 Z and y f = 7o + 5y;, using © and ©, we have 
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-2(ro + 8y r ) -1 
1 -2(70 + 57) 





\ ( x \ 
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yt 



where Cq = x\ + Vq = 1. It is clear that 



Hence, 



Q = lx t x t + ly t y t = -4(ft + 5y)(^ + y f 2 ) > -4(% + y)C t . 
q > e -4(%+r>. 



If f G [0, r"] where 



we have Q > C. 



rpll 



InC 



P 4(7o + 7) 



(82) 



(83) 



(84) 



(85) 



(86) 



I. Proof of Theorem U4\ 

Proof: For a single qubit system subject to dHJ, when H(t) = [1 + co(t)]I z + e x (t)I x + £ y (t)I y 
(\/ e x( t )+ £ y( t ) <£, £>0, \co(t)\ < ft) and co > 0), y = 7o + 5y (|5y| < 7), using ©, we have 



x t 



yt 



\ 



-4(7o + 57) -1 
1 -4(70 + ^7) 

-4(7o + 57) ) 



\ ( x \ 

Xt 



yt 



(87) 



where Pq = Xq + +Zq — 1 and R t — tr(p f 2 ) = xj+yj + zj. It is clear that 

Rt = -8(7o + 5%)P f > S(y + y)R t . 



(88) 



Hence, 

If t e [0, T d ] where 



we have P t > P. 

V. CONCLUSIONS 

Control design for quantum systems with uncertainties is an important task. This paper has 
proposed a sampled-data design approach for a single qubit with uncertainties. Both closed and 
Markovian open quantum systems are investigated, and uncertainties in the system Hamiltonian 
and uncertainties in the coupling strength of the system-environment interaction are analyzed. 
Several physically meaningful performance indices including fidelity, coherence and purity are 
used to define the required robustness and several sufficient conditions on the relationships be- 
tween related parameters in the control system are established to guarantee such robustness. The 
robust control law can be designed offline and then be used online on the single qubit system with 
uncertainties. Future work will include the extension of these sampled-data control approaches 
to other finite dimensional quantum systems and the development of practical applications of 
the proposed method. 
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